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Macrophytes are important sources of dissolved organic carbon (DOC) to littoral zones 
of lakes, but this DOC is believed to be mostly refractory to bacteria, leading to the 
hypothesis that bacterial metabolism is different in littoral and pelagic zones of a large 
subtropical shallow lake. We tested this hypothesis by three approaches: (I) dissolved inor- 
ganic carbon (DIC) accumulation in littoral and pelagic water; (II) O2 consumption estimate 
for a cloud of points (n = 47) covering the entire lake; (III) measurement of O2 consump- 
tion and CO2 accumulation in dark bottles, PCO2 in the water, lake-atmosphere fluxes 
of CO2 (fC02) and a large set of limnological variables at 19 sampling points (littoral and 
pelagic zones) during seven extensive campaigns. For the first two approaches, DIC and O2 
consumption were consistently lower in the littoral zone, and O2 consumption increased 
marginally with the distance to the nearest shore. For the third approach, we found in 
the littoral zone consistently lower DOC, total phosphorus (TP), and chlorophyll a, and a 
higher proportion of low-molecular-weight substances. Regression trees confirmed that 
high respiration (O2 consumption and CO2 production) was associated to lower concen- 
tration of low-molecular-weight substances, while PCO2 was associated to DOC andTR 
confirming that CO2 supersaturation occurs in an attempt to balance phosphorus defi- 
ciency of macrophyte substrates. Littoral zone fC02 showed a tendency to be a CO2 sink, 
whereas the pelagic zone showed a tendency to act as C0 2 source to the atmosphere. 
The high proportion of low-molecular-weight, unreactive substances, together with lower 
DOC and TP may impose lower rates of respiration in littoral zones. This effect of perennial 
stands of macrophytes may therefore have important, but not yet quantified implications 
for the global carbon metabolism of these lakes, but other issues still need to be carefully 
addressed before rejecting the general belief that macrophytes are always beneficial to 
bacteria. 
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INTRODUCTION 

Subtropical shallow lakes may present an important and differ- 
entiated set of conditions that make them ecologically distinct 
from other types of lakes. Because of the shallow mean depths and 
the benign climate conditions, macrophytes can colonize extensive 
areas and grow continuously throughout the year. These plants are 
important sources of organic matter to the littoral zones; in many 
cases they contribute more carbon than do the algae (Wetzel, 1992; 
Lauster et al., 2006) and in some systems, they sustain most of the 
bacterial production (Stanley et al., 2003). 

However, macrophyte -derived carbon is believed to be mostly 
in the form of high-molecular-weight polymer-like compounds 
(Bracchini et al, 2006), poor in nitrogen and phosphorus contents 
(Hessen, 1992), and therefore refractory to bacterial consumption. 
Studies finding that the presence of macrophytes can be detri- 
mental to bacterioplankton are starting to accumulate. Wu et al. 
(2007) found lower bacterial diversity in submersed macrophyte - 
dominated areas than in areas with no macrophytes in the large 
Lake Taihu, China. Rooney and Kalff (2003) surveyed nine lakes 



with different percentages of macrophyte coverage, and found a 
significant decrease in the bacterial respiration rate with increas- 
ing macrophyte coverage. In southern Brazil, bacterial metabolism 
and biovolume were found to be lower in the littoral (covered with 
macrophytes) than in the pelagic zones of subtropical shallow lakes 
(They et al., 2010). If bacterial metabolism can be lower in the pres- 
ence of macrophytes, there is a contradiction with the common 
belief that macrophyte -derived carbon always benefit bacteria, and 
it is reasonable to hypothesize that carbon cycling can be affected 
in the littoral zones of lakes that are extensively colonized by these 
plants. 

Estimates of bacterial respiration are essential to assess the 
balance of a lake's metabolism because they provide a direct assess- 
ment of the fate of carbon within the system, reflecting its trans- 
fer from the organic to the inorganic pool (Jahnke and Craven, 
1995). It is well established that plankton respiration increases 
with chlorophyll, phosphorus, and organic carbon concentrations 
in lakes (Pace and Prairie, 2005). The effect of primary produc- 
tion (in terms of chlorophyll a) on bacterial production has long 
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been recognized, even in systems dominated by terrestrial-carbon 
inputs (Kritzberg et al., 2005). Phosphorus can also play a role. 
High rates of respiration and CO2 supersaturation have been asso- 
ciated with high humic content of the carbon pool, because they 
are an attempt to balance phosphorus deficiency in these low- 
nutrient-content substrates (Hessen, 1992). Moreover, the ratio of 
phosphorus to labile DOC can control partitioning between bacte- 
rial respiration and production (Cimbleris and Kalff, 1998; Smith 
and Prairie, 2004; del Giorgio and Newell, 2012). DOC concentra- 
tion, in turn, is also positively associated with CO2 supersaturation 
in lakes (Prairie et al., 2002). 

In the present study, we tested the hypothesis that there is a dif- 
ference in the general metabolism of carbon in the littoral versus 
pelagic zones in the large, subtropical shallow Lake Mangueira, 
through: (a) an experiment to measure dissolved inorganic car- 
bon (DIC) accumulation in littoral and pelagic water; (b) mea- 
surement of respiration in a cloud of points covering littoral 
and pelagic zones; and (c) extensive temporal measurement of 
respiration, PCO2, CO2 fluxes, and limnological variables related 



to bacterial- carbon metabolism in a smaller cloud of points, also 
covering the whole lake. 

MATERIALS AND METHODS 
STUDY SITE 

Lake Mangueira (80,800 ha, southern Brazil, state of Rio Grande 
do Sul) is a large, shallow, freshwater subtropical coastal lake. It is 
surrounded by extensive belts of wetlands, located primarily in the 
northern and southern areas. The DIC accumulation experiment 
(to test question a above) was done in the southernmost part of 
the lake, and the two clouds of points (to test questions b and c 
above) covered the entire lake (Figure 1). 

DIC ACCUMULATION IN LITTORAL x PELAGIC WATER 

Samples of pelagic and littoral water were collected and incu- 
bated for 14 days for determination of the cumulative DIC, in 
April-May 2007. We retrieved vials at days 0, 1, 3, 5, 7, 10, and 
14, in triplicates (n = 42). The incubation vials (sterile 40 mL 
clear borosilicate vials, open top with silicone/PTFE septum and 
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FIGURE 1 | Locations of sampling points for respiration rates 
sampled in March 2009 (main map), and extensive sampling 
campaigns of respiration rates and limnological variables carried 



out in 2010-2011 (inset) in Lake Mangueira, a large shallow lake in 
subtropical southern Brazil. Maps are georeferenced in UTM 
coordinate system. 
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teflon-lined screw caps) were carefully rilled with no headspace 
(Farjalla et al., 2001) and incubated at room temperature (~20°C) 
in the dark. DIC was analyzed in a Total Organic Carbon (TOC) 
Analyzer (Shimadzu VCPH). Bulk respiration rates were consid- 
ered an approximation of bacterial respiration, since several pilot 
experiments with different water fractions showed no significant 
differences between the bulk and bacterial fractions after filtra- 
tion in MN 640d Macherey-Nagel paper filters (2.0-4.0 |xm mean 
retention size). This filtration has been extensively tested and 
successfully excludes ciliates, metazoans, and most flagellates. 

RESPIRATION IN LITTORAL x PELAGIC ZONES 

A large set of measurements of respiration rates, covering the entire 
Lake Mangueira, was taken in March 2009. The water was collected 
in the littoral and pelagic zones (littoral n = 2l; pelagic n = 26) and 
the water was incubated in dark bottles, for 6 days in the labora- 
tory (Wetzel and Likens, 2000). Dissolved oxygen was measured 
by a TOGA/SRI gas chromatograph (purge trap injector for water 
samples and TCD/HID detectors). Oxygen measurements were 
converted to carbon, assuming a molar conversion factor of 1.0 
(del Giorgio et al., 1997). Results were expressed in |xgC L _1 h _1 , 
assuming a constant respiratory rate. 

Nineteen points in the littoral and pelagic zones (littoral n = l2; 
pelagic n = 7) covering the entire lake were collected for sev- 
eral limnological variables during seven campaigns: pH, Abs250, 
Abs365, Abs250:365, total phosphorus (TP), total nitrogen (TN), 
chlorophyll a, DOC, DIC, O2 consumption, CO2 accumulation, 
PCO2 in the water, and CO2 fluxes (/CO2). With the exceptions of 
PCO2 and CO2 fluxes, all variables were collected in May, August, 
and November 2010 and March, May, August, and November 2011. 
PCO2 and /CO2 were estimated in February, May, August, and 
November 2010 and March, June, and August 201 1. 

Surface water was taken using a horizontal sampling bottle. pH 
was measured with a potentiometer (Wetzel and Likens, 2000). 
Abs250 and Abs365 refer to absorbance at 250 and 365 nm respec- 
tively, in 1 cm-optical path length quartz cuvettes. The former is 
proportional to the low- molecular- weight compounds, while the 
latter is proportional to the high-molecular-weight compounds 
present; and Abs250:365 is the ratio Abs250:Abs365 (Strome and 
Miller, 1978; Lindell et al, 1995). Samples for TN and TP were 
frozen in 1 L polyethylene bottles, and these nutrients were mea- 
sured by colorimetry (Mackereth et al., 1989). Chlorophyll a 
was quantified through cold ethanol extraction (Jespersen and 
Christoffersen, 1987). DOC and DIC samples were collected in 
30 mL precombusted (450°C for 1 h) amber glass bottles. DOC 
samples received a few drops of i^PO^ 3 , and DIC and DOC were 
analyzed in the same TOC analyzer (Shimadzu VCPH). 

Oxygen consumption and CO2 accumulation were determined 
by incubation in glass serum bottles. The bottles were kept inside 
a thermal box (approximately 23°C) and taken to the laboratory, 
totaling 120 h of incubation. O2 was measured in a TOGA/SRI 
gas chromatograph, and CO2 indirectly from measurements of 
DIC and pH after correction of temperature and ionic strength 
(Stumm and Morgan, 1996). The measurements were made at the 
beginning and at the end of incubation. The ionic strength was 
estimated from the electrical conductivity of the water measured 
in the field, using the following conversion formula (Snoeyink and 



Jenkins, 1980): 

IS ~ 1.6xl0 -5 x Conductivity 

where IS is the ionic strength and the water conductivity is 
measured in |xS cm -1 . 

The pCC>2 in the water was calculated based on the Law of 
Henry, with the following equation: 

where ao corresponds to the fraction of DIC as CO2 and Kh is the 
dissolution constant of CO2 expressed in moles L _1 atm -1 cor- 
rected for temperature in kelvin ( 7\) after the equation of Weiss 
(1974): 

In K H (C0 2 ) = -58.0931 + 90.5069 (^!jr) + 22 - 294 ln (^) 

The flux of CO2 (/CO2, mmolesCm -2 d _1 ) between atmos- 
phere and lake was estimated with the following equation (Mac- 
Intyre et al, 1995): 

fC0 2 = k x (pCOf m - pCO?) 

where pCOf m is the partial pressure of CO2 in the atmosphere (for 
380 ppm), pCO™ is the partial pressure of CO2 estimated for the 
surface layer of the lake, and k x is the coefficient of mass transfer 
(cm h _1 ), given by the following empirical relationship: 

kx = hoo {^) 

where x is equal to 0.66 for winds <3.0 m s _1 and equal to 0.5 for 
winds >3.0ms _1 , S c is the number of Schmidt for CO2, which 
is dependent on temperature (°C), according to the following 
relationship (Wanninkhof, 1992): 

S c (C0 2 ) = 1911.1 - 118.11T + 3.4527T 2 - 0.04132T 3 

where k^oo is estimated from wind velocity (Cole and Caraco, 
1998): 

^600 = 2.07 + (0.2 1 5 LJ/q 7 ) 

where 17" 10 is the wind velocity (m s _1 ). 

Data from the Santa Vitoria do Palmar meteorological station 
(www.inmet.gov.br), located 7.0 km from the lake, were used to 
obtain air temperature and wind velocity. 

STATISTICAL TREATMENT 

Differences between DIC accumulation in dark (littoral) and clear 
(pelagic) water samples and among different time periods for each 
type of water were tested with one-way non-parametric ANOVA 
(NPMANOVA in the univariate mode) test (Anderson, 2001). 
The dissimilarity measure used was Euclidean distance, with 9999 
permutations. Software: PAST 2.14 (Hammer et al., 2001). 
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FIGURE 2 | Respiration rates in littoral (n = 21) and pelagic (n = 26) 
zones (A) and regression of respiration rates against the distance to 
the nearest shore (B) in Lake Mangueira, sampled in March 2009. 



For the respiration data set, the mean rates of respiration in 
the littoral (n = 21) and pelagic zones (n = 26) were tested for dif- 
ference with ANOVA after natural-log transformation (correction 
for normality). The association between rates of respiration (also 
transformed) of the cloud of points and the distances from the 
nearest margin was tested by linear regression (Software R 2.15.0, 
R Development Core Team, 2012). 

The data from the O2 consumption, CO2 accumulation, pC0 2 , 
CO2 fluxes, and limnological variables obtained in the seven sam- 
pling campaigns were z-score standardized within each month and 
split into two categories, based on the distance from the nearest 
shore: littoral (< 1000 m) and pelagic (> 1000 m). This threshold 
was based on the mean distance of the samples from the shore. Dif- 
ferences between the littoral and pelagic zones were tested for each 
variable individually (pH, Abs250, Abs365, Abs250:365, TP, TN, 
chl a, DOC, 0 2 , C0 2 , jC0 2 ) by one-way NPMANOVA (Ander- 
son, 2001). We used the software PAST v. 2.14 (Hammer et al., 
2001) with Euclidean distance and 9999 permutations. A posteri- 
ori contrast tests used Bonferroni correction (n x p, where n is the 
number of comparisons and p is the p value obtained by permu- 
tation) implemented in the software PAST 2.14 (Hammer et al., 
2001). 

The explanatory power of abiotic variables (pH, Abs250, 
Abs365, Abs250:365, TP, TN, chl a, and DOC) on 0 2 con- 
sumption, CO2 accumulation, pC0 2 , and fC0 2 was assessed by 
binary recursive partitioning (regression tree) with the help of 
the R 2.15.0 (R Development Core Team, 2012) package {rpart} 
(Therneau et al., 2012) and following Crawley (2007). For O2 
consumption and CO2 accumulation all months were included 
(n = 133); for CO2 and /CO2, given the unbalance in sam- 
pling months, February (2010) and June (2011) were excluded 
(n = 95). 

RESULTS 

DIC ACCUMULATION EXPERIMENT 

The NPMANOVA test indicated that DIC was higher in the 
pelagic zone (23.24 ± 1.54 mgCL -1 ) than in the littoral zone 
(16.82 ± 1.72 mg C L~ l ) y pseudo-F = 163.0, p = 0.0001. There was 
no change in DIC with time: Littoral zone: pseudo-F = 0.957, 
p = 0.524; pelagic zone: pseudo-F = 0.685, p = 0.777. 

RESPIRATION CLOUDS OF POINTS 

The respiration data collected from the entire lake revealed that 
respiration rates were significantly lower in the littoral zone: 
F(l, 45) = 10.778, p < 0.002 (Figure 2A). The regression between 
respiration rates and distance from the nearest shore was margin- 
ally significant: ln(Respiration) = —2.648 + 0.559 [ln(Distance)], 
Adjusted R 2 = 0.046, p = 0.079 (Figure 2B). 

The NPMANOVA with all months revealed a significantly 
higher Abs250:365 in the littoral zone and a significantly higher 
TP and chl a in the pelagic zone. DOC was marginally significantly 
higher in the pelagic zone (Tables 1 and 2). 

The NPMANOVA by each month separately revealed that TP 
and chl a were frequently higher in the pelagic zones (May, August 
2010 and March, August, and November 201 1). Less frequent, but 
also consistent was higher Abs250 and Abs250:365 in littoral zones 
(May 2010 and May and August 2011). Punctually, littoral zones 



also showed higher pH (May 2010) and DIC (August 201 1), while 
the pelagic zone showed higher CO2 concentration (November 
2010; Tables 1 and 2). 

The regression trees revealed that the Abs250 was the major 
environmental variable influencing O2 consumption, CO2 pro- 
duction, and pC0 2 , while DOC influenced most the fC0 2 
(Figures 3A,B). 

For O2 consumption, the highest mean occurred at 
Abs250 < 0.064; at Abs250 > 0.064, there was also influence 
of Abs365, chl a, DOC, and pH, with high O2 consump- 
tion occurring at Abs365> 0.017 and DOC > 2.54 mgL -1 . At 
lower concentrations of DOC (<2.54mgL _1 ), O2 consump- 
tion was also high when the pH was >8.2. Lowest O2 con- 
sumption occurred at Abs365 < 0.017 and chl a concentrations 
>4.02|xgL _1 (Figure 3A). 
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FIGURE 3 | Explanatory power of environmental variables assessed by 
regression tree on 0 2 consumption (mg L _1 ) (A) and C0 2 (mg L 1 ) 
production (B) with the extensive campaigns data (/? = 133). Means in 
terminal leafs. Notice that the partitioning criteria above nodes applies to 
the left branches (< or >), while the complement (< or >) applies to the 
right branches. 



Table 2 | Non-parametric MANOVA (NPMANOVA) analysis on 
individual environmental variables for all campaigns and by month 
for the cloud of points sampled in lake Mangueria from 2010 to 2011. 



Year Month Variable NPMANOVA 







uutput 


ttrect 


2010-2011 All 


Abs250:365 


F* = 5.20, p = 0.023 


L>P 




TP 


F* = 4.02, p = 0.044 


P>L 




Chi a 


F* = 12.44, p< 0.001 


P>L 




DOC 


F* = 3.11, p = 0.083* 


P>L 


2010 May 


pH 


F* = 3.95, p = 0.066* 


L>P 




Abs250:365 


F* = 7.02, p = 0.0172 


L>P 




TP 


F* = 9.24, p = 0.010 


P>L 




Chi a 


F*=4.15, p = 0.062* 


P>L 


August 


Chi a 


F* = 4.78, p = 0.042 


P>L 


November 


C0 2 


F* = 3.14, p = 0.093* 


P>L 


2011 March 


TP 


F* = 2.91, p = 0.098* 


P>L 




Chi a 


F* = 3.06, p = 0.090* 


P>L 


May 


Abs250 


F* = 3.67, p = 0.042 


L>P 




Abs365 


F* = 3.29, p = 0.044 


L>P 


August 


Abs250 


F* = 4.03, p = 0.045 


L>P 




Abs250:365 


F* = 4.41, p = 0.013 


L>P 




Chi a 


F* = 3.97, p = 0.0597* 


P>L 


November 


TP 


F* = 7.14, p = 0.015 


L>P 


Groups coded according 
(L, < 1000 m) x Pelagic (P, 


to the distance from the nearest shore: littoral 
x > 1000m). In bold: p<0.05; * 0.05<p<0. 1). F\ 



pseudo-F. 

CO2 production was also high at Abs250 < 0.064, but 
reached higher means at Abs250 > 0.064), Abs250:365 < 7.35, 
DOC < 2.51 mgL -1 , pH > 7.67, and high chl a concentra- 
tions (between 2.59 and 3.07 |xgL -1 ). At DOC concentra- 
tions <2.51mgL _1 and low pH < 7.67, the mean was also 
high. The low values, however, occurred at Abs250 > 0.064, and 
Abs250:365 > 7.35 (Figure 3B). 

pC0 2 showed high values at Abs250 > 0.104. At lower Abs250 
(< 0.104), there was also influence of DOC, TP, pH, and 
Abs250, with lowest values occurring at DOC between 0.58 and 
3.95 mg L" 1 , TP > 0.029 mg L" 1 , and pH > 7.70 (Figure 4A). 

JCO2 showed a strong relationship with DOC concentra- 
tions. At DOC >2.54mgL _1 , and Abs250 > 0.0853 it was found 
the highest positive value of /CO2. A negative mean was found 
with values of DOC between 1.48 and 2.54 mgL -1 and at 
Abs250 > 0.0846 (Figure 4B). 

DISCUSSION 

Our results showed that respiration can be lower in the littoral 
zone, based on the consistently lower DIC in water collected from 
this zone in incubation assays. In the first whole-lake survey, we 
found significantly lower respiration rates in these zones and a 
marginal (p = 0.079) increase of respiration rate with distance to 
the nearest shore. In the extensive sampling, this was not consis- 
tently found throughout all months. Marginally (p = 0.093) lower 
CO2 accumulation was detectable in at least 1 month (Novem- 
ber 2010) in the littoral zones. However, significantly higher DIC 



occurred in the littoral zones in August 201 1. These differences are 
possibly associated with differences in environmental conditions 
between these zones, as shown by consistently higher chlorophyll a 
and TP, marginally higher DOC (p = 0.083) and significantly lower 
Abs250:365 in the pelagic zone, as shown by NPMANOVA analy- 
ses on the entire data set (all months). The regression tree analyses 
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FIGURE 4 | Explanatory power of environmental variables assessed by 
regression tree on pC0 2 (|iatm) (A) and fC0 2 (mmolCm 2 d _1 ) (B) with 
extensive campaigns data (n = 95). Means in terminal leafs. Notice that 
the partitioning criteria above nodes applies to the left branches (< or >), 
while the complement (< or >) applies to the right branches. 



confirmed that higher O2 consumption and CO2 production are 
associated to low Abs250, high DOC, and high Abs365. PCO2 and 
/CO2, at their turn, were higher at high Abs250, DOC, and TP. 

Our results showed contrasting patterns. The DIC accumula- 
tion experiment and the extensive sampling of 2009 supported 
our hypothesis that respiration is lower in the littoral zone than in 
the pelagic zone. However, apart from the marginally higher CO2 
accumulation in November 2010 in the pelagic zone, the extensive 



sampling of 2010-2011 did not show such striking evidence in 
support of the hypothesis. This could be due to the smaller num- 
ber of points sampled, but also the result of intensive mixing of 
Lake Mangueira by wind (Cardoso et al, 2012), which makes it 
difficult to detect abrupt differences in respiration in the field, as 
indicated by several marginal relationships found. This could have 
been crucial because of the smaller set of sampling points. 

One hypothesis for the lack of detection of DIC accumulation is 
that in areas with an extensive coverage of macrophytes, like south- 
ern Mangueira, the macrophyte influence might reach pelagic 
zones as well. We unfortunately do not have enough data gathered 
yet to envisage how, why, and when this occurs. We nevertheless 
did observe that DIC was on average higher in the pelagic zone, 
likely due to very low accumulation rates (not statistically signifi- 
cant in < 14 days). When we extend our perspective (spatially and 
temporally) to the whole lake and for different years and seasons, 
we see a different picture, where we do detect respiration (even 
though still low) enough to put in evidence differences between 
littoral and pelagic zones. Another important point to take into 
consideration is that lake Mangueira is of marine origin (former 
closed estuary) with significant accumulation of sea shell banks. 
This explains the high background levels of DIC in the lake, and 
could have masked measurements of bacterial respiration based on 
DIC. This also suggests that it is not the best method for this pur- 
pose, and may have been the cause for the disagreement between 
the patterns of DIC and O2/CO2 respiration rates observed. 

As stated, in Lake Mangueira, rates of respiration are 
extremely low, and it is striking that it is impossible to 
detect respiration in dark bottles over a period shorter than 
5 days. The mean respiration rates in the littoral (1.34mmol 
02m _3 h _1 ) and pelagic zones (1.43mmol 02m _3 h _1 ) are 
closer to the lower range presented in the extensive compilation 
of plankton respiration rates (0.029 to 6.73 mmol 02 m -3 h _1 ) 
by Pace and Prairie (2005). This is also confirmed in a 
carbon-basis comparison, as the rates of respiration in this 
study (Littoral zone: Mean: 1.469 |xgC L _1 h _1 , range: 0.008- 
6.241 ixgCL"^ -1 ; Pelagic zone: Mean: 6.540 |xgCL _1 h _1 , 
range: 0.250-23.718 |xgCL _1 h _1 ) are also close to the lower range 
reported in other systems: Delaware Bay and salt marshes (range: 
2.53-13.14 ixgCL^h- 1 ; del Giorgio and Newell, 2012) and 20 
Quebec lakes (range: 0.168-2.138 |xgCL _1 h" 1 ). 

The differences in respiration in the littoral and pelagic zones 
of Lake Mangueira may reflect differences in the responses of 
bacterial metabolism to the following drivers: (I) quantity and 
bioavailability of DOC; and (II) phosphorus limitation. Our 
results supported both hypotheses, as indicated by the lower chl a, 
DOC, and TP, and the higher proportion of low-molecular- weight 
compounds (in terms of Abs250 and Abs250:365) in the littoral 
zone found in many cases. 

The generally low DOC bioavailability in the form of humic 
substances (Minister and Chrost, 1990) indicates that a large pro- 
portion of dissolved carbon in Lake Mangueira is refractory to 
bacterial consumption. When compared to the values reported by 
Amado et al. (2006) for an Amazonian river and a stream, the con- 
tent of humic substances in Lake Mangueira is 6 times higher. The 
humic content is at least 4 times higher than for values reported 
for the mesohumic Lake Sjattesjon (Lindell et al, 1995), or even as 
much as 12-17 times higher than reported for the Gutierrez River 
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(all comparisons on the basis of the Abs250:365 nm ratio; Perez 
etal.,2003). 

DOC is generally high in densely vegetated areas (Wetzel, 1992; 
Reitner et al, 1999); however, we found higher DOC in the pelagic 
zone, along with chl a, suggesting the importance of algal car- 
bon to the system. This is in accordance with the long-recognized 
bacterial dependence on algal production (Cole, 1982), particu- 
larly because the more labile algal DOC sustains higher bacterial 
growth efficiency (Kritzberg et al., 2005). This was evident by 
the greater phytoplankton biomass concentrated (higher chloro- 
phyll a) in the pelagic zone at the end of winter (August 2010), 
which was respired by bacteria with the rising temperatures in 
spring (highest CO2 accumulation found in November 2010; 
Table 1). 

The higher concentration of low-molecular- weight substrates 
in the littoral zone can also imply less availability of substrates 
(Amon and Benner, 1996). Contrariwise, these substrates are sup- 
posedly more easily taken up by bacteria, but they may be more 
refractory to consumption. Tranvik (1990) found higher bac- 
terial production per unit of carbon in high-molecular-weight 
DOC than in low- molecular- weight DOC in lakes with different 
humic content in Sweden. In another study, Amon and Benner 
(1996) found higher bacterial production and respiration in high- 
molecular- weight (HMW, > 1 KDa) than low-molecular- weight 
(LMW, < 1 KDa) DOC. Based on these results, the authors pro- 
posed the size-reactivity continuum model, which predicts that 
the major path of degradation goes from large, highly reactive to 
small, highly recalcitrant molecules. If this hypothesis is broadly 
applicable to many ecosystems, it suggests that in Lake Mangueira 
the compounds lixiviated by the macrophytes may undergo some 
degradation, but accumulate as dissolved, unreactive LMW com- 
pounds with very low degradation rates in the littoral zones, in 
agreement with the higher Abs250:365 and lack of change in the 
DIC with time in the incubation experiment. Our results strongly 
support this hypothesis, as higher O2 consumption, and CO2 pro- 
duction were found for low Abs250 and high Abs365 or conversely 
lower O2 consumption, and CO2 production were found when the 
ratio Abs250:365 was too high. 

The lower TP found in the littoral zone can be attributed to the 
generally low watershed load and also the great competitive advan- 
tage of macrophytes in taking up this nutrient. As noted by Vidal 
et al. (2011), not only carbon, but the carbon: phosphorus ratio 
controls bacterial production and respiration in lakes, and even a 
situation of labile DOC accumulation can occur under phospho- 
rus limitation (Vadstein et al, 2003). In fact, increased bacterial 
density in the presence of macrophytes after phosphorus amend- 
ment has been described (Huss and Wehr, 2004; Joniak et al., 2007; 
Morozova et al, 2011). An interesting result was the lower TP in 
the pelagic zone in November 201 1, which could be due to deple- 
tion of this nutrient in the pelagic zone during summer (Table 1); 
however, no similar depletion occurred in 2010, and cannot be 
assumed as a pattern based on our data. The regression trees indi- 
cated an effect of TP when DOC is high; low TP is associated 
to low CO2 production, while the opposite situation also applies, 
supporting the hypothesis that low respiration can be associated 
to low TP. This is expected to generate lower pCC>2 in littoral zones 
where TP is lower; however a difference in pCC>2 between littoral 



and pelagic zones was not detected. The regression tree revealed 
that PCO2 was low when TP was high in samples of high DOC. 
This is in agreement with low phosphorus leading to CO2 super- 
saturation in order to balance this nutrient deficiency (Hessen, 
1992), and suggests that this also occurs in lake Mangueira. As 
CO2 concentration is only part of PCO2, other factors affecting 
PCO2 equilibrium need to be taken into account, and this may be 
the cause for the apparent contradiction between CO2 production 
and PCO2 predicted from TP. 

Levels of pCOj were in general lower than in other lakes 
with similar or even lower TP content: mean around 660 |xatm 
(range: 130-1010 |xatm) for a range of Quebec lakes with vari- 
able trophic content; mean around 1 036 |xatm (range: 1-20 
249 |xatm) for a comprehensive compilation of data for 1835 lakes 
distributed worldwide (Cole et al, 1994). These levels are far above 
those found in this study. This could be explained by the gener- 
ally low DOC concentration in Lake Mangueira (mean around 
2.2 mgL -1 ), because it has been demonstrated that lakes are net 
heterotrophic (supersaturated with CO2) above concentrations of 
4-6 mgL -1 DOC (Prairie et al, 2002). This contrasts with the 
generally accepted view that most aquatic systems are net het- 
erotrophic and function as net sources of CO2 to the atmosphere 
(Cole et al, 1994; Duarte and Prairie, 2005), which makes Lake 
Mangueira (particularly the littoral zone) an exception to this 
general rule. 

The analysis of the frequency of occurrence (histogram) of 
classes of /CO2 between the lake surface and atmosphere in the 
littoral zone showed a higher frequency of positive (net influx of 
CO2 to the lake), while the pelagic zone showed a higher frequency 
of negative values (net efflux of CO2 from the lake) of /CO2, even 
though there was no difference between the means in the littoral 
and pelagic zones. A closer inspection of the fluxes showed that 
they are very close to zero, and the lake functions as a sink of 
CO2, especially during winter and spring (May to November), 
when the/C02 becomes close to zero or even positive (Figure 5). 
The regression tree showed that high, positive /CO2 was asso- 
ciated to higher DOC and Abs250, contrary to the expectation 
that higher DOC is associated to higher CO2 supersaturation and 
hence net efflux of CO2 from the lake; negative /CO2 was associ- 
ated to intermediate values of DOC and high Abs250. However, 
more direct evidence needs to be gathered in order to draw more 
solid associations. 

Besides main effects of DOC and nutrients, other possible 
explanations for low respiration in the littoral zone in Lake 
Mangueira have also been hypothesized: formation of hydrogen 
peroxide by UV radiation, release of allelochemicals by macro- 
phytes, and limitation by micronutrients, e.g., iron and silicate. 
UV exposure of macrophytes products of degradation in Lake 
Mangueira has been demonstrated to decrease bacterial produc- 
tion, possibly by formation of H2O2 (They et al., submitted). The 
release of allelopathic substances by (mainly submersed) macro- 
phytes is widely known (Gross et al., 2007; Mulderij et al, 2007), 
and has been considered as a possible explanation for lower bac- 
terial diversity (Wu et al., 2007) and metabolism (They et al., 
2010) in areas extensively colonized by macrophytes. We have data 
on Fe and silicate from the extensive campaigns that were not 
included in the results because they are out of objectives of this 
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FIGURE 5 | Histogram of C0 2 fluxes between atmosphere and lake 
surface in littoral (n = 12) and pelagic (n = 7) zones in Lake Mangueira, 
sampled in 2010-2011. 



study. Based on these data, there is no evidence of iron limitation, 
since the values found in the littoral zones (0.325 ± 0.253 mgL" 1 ) 
and pelagic zones (0.349 =b 0.248 mg L _1 ) are far above concentra- 
tions reported to be limiting in other systems, like, e.g., open ocean 
(<5.6x 10 -6 mgL -1 ; Oliver etal, 2004) or lake Erie (1.7 x 10 -4 - 
1.09 x 10~ 2 mgL -1 ; North etal., 2007). Silicate concentrations are 
far above those reported to be limiting (2 |xm) for diatoms (Egge 
and Aksnes, 1992) in littoral (100.92 ± 23.99 |xm) and pelagic 
zones (100.04 zb 23.83 |xm), and thus there no evidence of possi- 
ble impacts on bacterial supply of organic substrates derived from 
silicate limitation of phytoplankton. 

One important and final consideration is that we found support 
for our hypotheses by different, independent methods: one-time 
DIC accumulation, extensive one-time sampling and extensive 
many-times sampling by O2 consumption, and CO2 production 
measurements. Each method has its own limitations, as we 



believed was crucial in the DIC experiment; had we employed 
a longer incubation period, we could have seen some detectable 
accumulation. Differences between littoral and pelagic zones has 
long been recognized, but what is generally believed is that macro - 
phyte presence/carbon is always good to bacteria, a paradigm our 
data do not support. 

The littoral zone of Lake Mangueira (especially in the north 
and south) is extensively colonized by emergent and submersed 
macrophytes, respectively, and these plants are expected to con- 
tribute large amounts of organic carbon to the system. However, 
our results showed that respiration in this zone can be lower 
than in the pelagic zone, at least during a part of the year. This 
may be the result of differences in DOC quality, mainly derived 
from macrophytes in the littoral zone and phytoplankton in the 
pelagic zone. The disappearance of these differences may be due 
to seasonality or to masking by mixing in this shallow, wind- 
dominated lake. The finding of lower respiration rates in the 
littoral zone means that DOC remains in the system, mostly in a 
low-molecular- weight, unreactive form. Hence, the general belief 
that macrophyte- derived carbon is always beneficial to bacteria 
is not supported. The littoral zone, therefore, shows a greater 
tendency to be a CO2 sink, compared to the pelagic zone. If 
lower respiration in littoral zones is a common feature of sub- 
tropical shallow lakes dominated by macrophytes, there may be 
important and still unrecognized implications for their global 
carbon metabolism. However important these implications may 
seem, important issues (e.g., organic carbon molecular size and 
quality spectra, effect of taxonomic structure of macrophytes, bac- 
terial taxonomic and functional diversity, the role of sediment 
bacteria, and more experimental and field evidence) need to be 
addressed for a better understanding of the many ways the interac- 
tion between macrophytes and bacteria can impact global carbon 
metabolism of lakes. 
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